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Abstract 

We present a first-principles study of native defects in NaAlH4. Our analysis indicates that the 
structure and energetics of these defects can be interpreted in terms of elementary building blocks, 
which include V^J H4 , V^ a , V^, H~, and (H^j. We also calculate migration barriers for several key 
defects, in order to compare enthalpies of diffusion to experimentally measured activation energies 
of desorption. From this, we estimate activation energies for diffusion of defects and defect pairs. 
We suggest that V^J H4 and H7, or V^7 a and V^, may be responsible for diffusion necessary for 
desorption. We discuss the possible role of V^~-H~ complex formation. The values we find are in 
the range of activation energies reported for catalyzed desorption. 
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I. INTRODUCTION 



NaAlH4 is an interesting hydrogen storage material. While its theoretical hydrogen ca- 
pacity by weight (5.6 %) is not sufficient for automotive applications, it may be useful in 
other applications. More importantly, as one of the most widely studied hydrogen storage 
materials, it serves as a prototype for fundamental investigations of kinetics. The struc- 
ture is shown in Fig. [TJ It can most easily be understood as an NbP ordering of Na + and 
(AIH4) - tetrahedra. Under practical conditions, NaAlH-4 undergoes a two-step reaction to 
release hydrogen: 

NaAlH 4 -Na 3 AlH 6 + -Al + H 2 (1) 
3 3 

^Na 3 AlH 6 -> NaH + ±A1 + ^H 2 . (2) 

Reversible absorption and desorption at reasonable temperatures was first accomplished 
by Bogdanovic and Schwickardi in 1997 by adding a few percent titanium. [jj The mechanism 
of this kinetic improvement has remained controversial. Recent first-princples calculations 
suggest that titanium may play a role as an electronically active impurity, promoting the 
diffusion of hydrogen. |2j A number of kinetic experiments have been performed on the 
above reactions, 3, 4, and these studies are in reasonable agreement with regards 

to the reported activation energy of desorption. Desorption of hydrogen and decomposi- 
tion of NaAlH-4 requires not only mass transport of hydrogen but also of aluminum and/or 
sodium. This process is likely to be mediated by native defects. Lohstroh and Fichtner 
suggested that desorption from NaAlH-4 is rate-limited by diffusion. 7J Diffusion of native 
defects was recently approached by Gunaydin et al. through first-principles methods. 9| 
Here we devote further study to this topic. 

In this paper we investigate structure and stability of native defects in NaAlH 4 based on 
first-principles density-functional theory. For relevant defects, migration enthalpies are also 
calculated. These allow us to estimate diffusion activation energies for the various defects 
that may be responsible for mass transport. We find that most of the relevant defects exist 
in charge states other than neutral, and that consideration of these charge states is essential 
for a proper description of migration and kinetics. Section [TT] describes the computational 
approach. In Sec.fllllwe report our results. Section [TV] contains a discussion and comparison 
with experiment. 
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II. METHODS 



We use the Vienna Ab-initio Simulation Package (VASP) |10l . I 111 . Il2j | to perform den 



sity functional theory (DFT) calculations within the Perdew-Berke-Erzerhoff (PBE) [13 1 
generalized gradient approximation (GGA). The plane-wave cutoff is 500 eV, and projector- 
augmented wave potentials [ijj] are used. In the case of sodium, 2p orbitals are included 
in the valence description of the atom. Our calculated theoretical lattice parameters for 
NaAlH4 are a = 5.01 A and c = 11.12 A, within 0.7% of the experimental values, [ijj] 

We calculate point defects in a supercell containing 96 atoms. The supercell dimen- 
sions are kept at the theoretical bulk lattice parameters, but of course the atoms within 
the cell are fully re 
aluminum hydride. 



axed. This method has been used previously in the study of sodium 
2], The formation energy of a defect is a key quantity, determining 



its concentration in the lattice through the relation: 17J] 



c(X) = iV sitcs iV config exp[E f (X)/kT], (3) 

where E^(X) is the formation energy of defect X, N s it es is the number of lattice sites per unit 
volume on which the defect can be incorporated, and A^ on fig is the number of configurations 
per site in which the defect can be formed. 

NaAlH4 is an insulator with a wide band gap - it is expected that native defects exist in 
charge states other than neutral. There exists a clear prescription for calculating the forma- 
tion energies of charged defects in a manner that accurately describes the thermodynamic 



17( The formation energy of a charged defect X is 



reservoirs for atoms and for electrons, 
calculated as: 

^(X 9 ) = E tot (X 9 ) - £ tot (bulk) + qe F + mm- ( 4 ) 

i=0 

E tot (X <? ) is the total energy of a supercell containing the defect in charge state q. Similarly, 
E tot (bulk) is the total energy of a supercell without a defect. The last term in the expression 
ensures stochiometric balance, with representing the number of atoms exchanged with the 
reservoir with chemical potential /ij. Our chemical potentials are referenced to the standard 
state (i.e., bulk bcc Na, bulk fee Al, and H2 molecules at T=0). Presented separately 
from the atomic reservoir chemical potentials is ep, the chemical potential for electrons or 
Fermi energy. We keep with convention and reference the Fermi energy to the valence-band 
maximum of the bulk material. In this study, a correction to the formation energy of a 



charged defect is made through averages of the electrostatic potential in regions far away 
from the defect, [l^] 

Migration energies were calculated by using an implementation of the Nudged- Elastic 
Band (NEB) method [3] within VASP. For charged defects, in principle the aforementioned 
correction related to the electrostatic potential alignment could vary along the diffusion 
path; we have checked that for singly charged defects this correction term is small. In the 
case of diffusion of V^ a , inclusion of this energetic correction caused a decrease of 0.06 eV in 
the energy of the saddle point relative to the minima. 

The enthalpy of diffusion is the sum of formation and migration enthalpies. These en- 
thalpies do not include any pressure term from gas-phase H 2 - to first approximation this 
term should drop out of an activation energy measurement. In the following, we therefore 
focus on activation energies for diffusion that we take to be the sum of calculated formation 
and migration energies. 



III. RESULTS 



A. Chemical potentials 

Our calculated formation energies are completely general, and can be applied to any con- 
dition described by a set of atomic chemical potentials. I.e., the atomic chemical potentials 
/ii in Eq.( H]) are variables that can describe different sets of experimental conditions. It is 
useful to consider various possible scenarios that lead to specific constraints. 

Equilibrium with NaAlH4 implies that 

/"Na + /iAi + 4/i H = AH f (NaAlH 4 ), (5) 

H 4 ; our calculated value for 



where Aiff(NaAlH 4 ) is the enthalpy of formation of NaA 

this quantity is —0.824 eV (experiment: —1.205 eV, Ref. Il9l ). Similarly, equilibrium with 
NasAlHg implies that 

3/iNa + /iAi + 6/i H = A#/(Na 3 AlH 6 ), (6) 
formation of Na-HAlH fi ; our calculated value is 



where Aifj(Na 3 AlH 6 ) is the enthalpy o 
-1.784 eV (experiment: -2.475 eV, Ref. 
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The difference between our calculated enthalpies of formation and experiment can be 
partly attributed to zero-point corrections to the vibrational energies. We do not con- 
sider such corrections here; their inclusion would not change our qualitative conclusions for 
formation energies, since typically significant cancellation occurs between terms in the defect 
and in the reservoirs, [l^] I.e., energy differences between comparable solids are typically 
well described. We note, for instance, that our calculated heat of reaction for Eq. is 0.23 



eV, and for Eq. © it is 0.20 eV; t 
0.38 eV and 0.24 eV, respectively. 



rese values are quite close to the experimental values of 



19[ 



For calculating defect concentrations, in principle one would have to include the contri- 
butions from finite temperature vibrational entropies. A rigorous treatment of this effect 
requires an evaluation of the vibrational spectrum for each of the defects. This would in- 
volve an extraordinary computational effort. Nonetheless, it would not lead to a deeper 
understanding of the effects discussed here, so we consider this task beyond the scope of the 
present work. Reasonable estimates indicate that such contributions are relatively small, 
mainly due to cancellation when comparing the vibration entropy of the defect system with 
the ideal host and reservoir. [itJ In addition, at the temperatures of interest in the present 
work (noting that NaAlH 4 is stable up to 200°C), the TS term will not make a significant 
contribution to the free energies of formation. 16] 

For purposes of presentation of formation energy results, it is convenient to choose a 
specific set of chemical potentials that are intended to be close to those relevant for dehydro- 
genation. We will choose the chemical potentials of Al, Na, and H by assuming equilibrium 
with NasAlHg, NaAlH-4, and Al. Using our calculated energies, this results in a value for 
/ih=— 0.12 eV, quite close to the Gibbs free energy of H 2 gas at 1 atm and 303 K, the equilib- 
rium temperature of NasAlUg, Al, NaAlH4, and H2 at 1 atm. This agreement supports 
our argument that this set of chemical potentials is representative of real conditions. 

While this particular set of chemical potentials presents a convenient and relevant set for 
presenting our formation energy results, it does not preclude us from examining situations 
that correspond to different choices of chemical potentials. The corresponding formation 
energies can easily be derived based on the values given in this paper and the general 
expression for formation energy. Examples will be given in Sec. IIVI 
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B. Hydrogen-related defects 



We start by presenting our first-principles results for h ydr ogen- related defects in NaAlH4 



in Fig. [21 These values have been calculated before, |2|, ll6j but we include them here for 



: comparison. Our present numbers are in good agreement with 



completeness and ease o 

the previous results; j^, 16j small differences arise from the use of a different exchange- 
correlation functional. In agreement with the previous calculations Q we find that the 
neutral defects, and H°, have higher formation energies than the charged defects for all 
Fermi-level positions. They are therefore not included in Fig. [2j 

H t ~ can be thought of as the addition of H~ to the system. The structure of H~ is 
shown in Fig. [3]^a). This defect can be viewed as an (A1 2 H 9 ) 3- unit; i.e., it is composed 
of two (AIH4) - units with an additional H~ located midway. The central hydrogen of this 
defect structure sits in a "bridge bond" -type arrangement between two aluminum sites at 
distances of 1.78 and 1.84 A. An isosurface of the charge density associated with the defect 
state, which occurs at 2.5 eV above the valence-band maximum, is included in Fig. [3](a). 
One might expect that, due to Coulomb interaction, H~ would prefer to sit next to a Na 
atom in NaAlH 4 . Indeed, we found a configuration in which Hr is located near two Na, with 
Na-H distances of 1.95 A. However, this configuration is metastable, being 0.6 eV higher in 
energy than the lowest energy configuration shown in Fig. [3](a). 

The migration path for this defect can be thought of as diffusion of the defect complex 
along the 4i screw axis of the system. One of the (AlhLi) - units to which the is attached 
slightly rotates and moves closer to an adjacent (AIH4) - until a new bridge bond forms, 
while simultaneously the bridge bond with the original (AIH4) - is broken. The resulting 
migration barrier of this defect is 0.16 eV. Such low values indicate a very high diffusivity 
of the point defect. Note that the migration barrier reported in the present work is slightly 
lower, by 0.06 eV, than the value reported in Ref. [J This small difference is attributed a 
slight difference in the migration paths. 

V^" can be viewed as the removal of H~ from the system. In response to this, the resulting 
undercoordinated Al site shifts towards an adjacent (Alh^) - tetrahedron. The defect can 
be thought of as a (A1 2 H 7 )~ complex [see Fig.JMb)]. The migration barrier for this species 



has previously been calculated to be 0.26 eV. 



2] 



Vfj can be thought of as the extraction of a proton from the system. Compared to 
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the large atomic rearrangements observed for the other defects, the formation of has 
a relatively small impact on the geometry [see Fig. [3](c)]. A plot of the highest occupied 
Kohn-Sham state (1.2 eV above the valence-band maximum) shows that an aluminum lone 
pair has replaced the missing H + . We have calculated the migration barrier for this defect by 
moving a H atom from a neighboring (A1H 4 )~ unit into the vacancy. The resulting migration 
barrier of this defect is 0.92 eV. This barrier is significantly higher than that of the other 
hydrogen-related defects. The reason is that the saddle-point configuration consists of a 
hydrogen atom located midway between two A1H 3 units. Such a configuration is favorable 
in the case of a positive charge state (see the case of but quite high in energy in the 
case of a negative charge state, due to the fact that an antibonding state resulting from the 
interaction between the Al atoms needs to be occupied in order to accommodate the charge. 

(H 2 )j is an interstitial molecule inside the hydride [see Fig. 12(d)] . Figure [2] shows that 
this defect is relatively high in energy. The calculated H-H bond length is 0.76 A, very close 
to that calculated for an isolated molecule (0.75 A). Migration of this defect has a calculated 
barrier of 0.25 eV. The bond length of the H2 dimer itself is very well preserved along the 
migration path. 

H^, finally, has a surprising structure [see Fig. [3](e)] that can be understood as a complex 
between and (H 2 )j, both discussed earlier [see Figs. Mb) and (d)]. Within this model, 
it is interesting to consider the binding energy of this pair. The formation energy of 
is lower than the sum of the formation energies of Vr and (H 2 ); by 0.71 eV. This binding 
energy is independent of any chemical potential for either atomic species or electrons. This 
relatively large binding energy can be attributed to the unfavorable sterics of the component 
(H 2 )j. When (H 2 )j is in the presence of a vacancy, the induced local strain in the lattice 
largely disappears. 

C. Aluminum-related defects 

Next we turn to native defects related to the Al site. The calculated formation energies 
are shown in Fig. HI 

V^ H4 corresponds to the removal of an entire (Alty - tetrahedron from the system. 
Figure [5]^a) shows that there is remarkably little structural relaxation around this defect. 
Figure H] shows that the formation energy of this defect is quite low, and we will see that it 
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plays an important role. NEB calculations for the migration path are presented in Fig. [61 
Migration consists of moving a neighboring tetrahedron into the vacancy. This tetrahedron 
remains fairly rigid along the migration path; see Fig. [6] for snapshots of the initial, saddle- 
point, and final configurations and Ref. 21 for an animation of the diffusion. This path 
results in a migration barrier of 0.46 eV. 

[Fig. [3(b)] can be understood as a complex of a V^ m and two (H 2 )j. The energy of 
V^J is lower than the sum of its components by 1.52 eV. As discussed earlier in the case of 
H~l , this can be explained on the basis of the unfavorable sterics of (H 2 )j in the absence of a 
vacancy. Indeed, we note that the 1.52 eV binding energy of V^, which contains two (H 2 )j, 
is very close to twice the binding energy (0.71 eV) calculated for H+, which contains one 
(H 2 )j species. 

V^f is a complex of V^ H4 and four rF~, as can be seen in Fig. [5](c). The binding energy of 
V^f, relative to the 5 component defects, is 1.46 eV. Even though the H7" component defects 
in V^f all have the same charge, Coulomb repulsion between them does not seem to play a 
major role. The fairly large binding energy of the V^f complex can again be attributed to 
sterics. 

V^j can be understood as a complex of V^j H , (H 2 )j, and two H,^. The binding energy is 
1.75 eV. The structure of is presented in Fig. [5](d). 

Vaih 3 [Fig. Ete)], finally, can be regarded as a complex of V^ lHi and H~. The binding 
energy, i.e. the difference in energy with respect to the sum of formation energies of these 
two constituents, is 0.37 eV. This binding energy is roughly one-quarter the binding energy 
of V^f, as one might expect - only one H~ is involved in Vaih 3 , whereas four H~ are used 
to construct V^f. This again implies that Coulomb repulsion between H~ does not play 
an important role, and that the binding energies that are gained when defects are brought 
together are primarily due to sterics. 

For completeness, we have also performed calculations for Al interstitials. We find, how- 
ever, that these defects all have formation energies that are significantly higher than the 
vacancy defects discussed above. The most favorable charge state is Al^ + , which has a for- 
mation energy of 4.18 eV at 6^=3.26 eV. This value is very high and therefore the defect is 
not included in Fig. [U 
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D. Sodium-related defects 



Formation energies for Na-related native defects are presented in Fig. [7J corresponds 
to the removal of a Na + ion from the system. As in the case of V^j H , there is remarkably 
little structural relaxation, as can be seen in Fig. [H](a). We calculated a migration energy of 
0.41 eV for this defect. 

VN a H is a complex of the defects V^ a and V^, as can be seen in Fig. [8]^b). The calculated 
binding energy is 0.14 eV. This relatively small binding energy is to be expected, since both 
component defects are vacancies and little energy can be gained due to strain relaxation. 

We have also investigated Na interstitials. Not surprisingly, Na^~ is most favorable; it has 
a formation energy of 1.44 eV at e^=3.26 eV and a migration barrier of 0.48 eV. Still, Na 
interstitials are less favorable than other defects. 



IV. DISCUSSION 

A. Comparison with previous calculations 

In Sec. IHIt we presented and compared formation energies of various native defects. 
In addition to the results presented in the figures, we also summarize key information for 
relevant defects in Table [H We noted that interstitial defects create local stress, and that this 
stress can be relieved by joining interstitials and vacancies. The binding energies between 
je explained by this sterics argument. 



such defects could 

Lodziana et al. 22j have reported first-principles calculations for native point defects in 
NaAlELi and LiBH4, using a methodology very similar to ours, though with a different choice 
of GGA exchange and correlation potential. Our own tests have indicated that the choice of 
GGA potential has only a minor impact on the results. Lodziana et al. also make a different 
choice of chemical potentials for presenting formation energies, invoking equilibrium with 
A1 2 H 6 , but again this should affect the calculated formation energies only by a few 0.1 
eV compared to our choice. Another difference is the use of supercell-size corrections for 
charged defects. We feel it is better not to apply uncontrolled approximations in an attempt 
to correct for such effects, since such corrections often tend to "overshoot" and make things 



worse rather than better 
results and those of Ref. 



Although there is broad qualitative agreement between our 
22j, a number of notable differences appear. For instance, the 



transition level between the + and — charge states of hydrogen vacancies and hydrogen 
interstitials occurs at much lower values in the band gap than it does in our calculations; the 
difference seems comparable to the energy of the valence-band maximum (highest occupied 
eigenstate) in the bulk ( —0.9 eV), which could indicate an error in the reference energy 
for the Fermi level. Lodziana et al. also find a significantly higher energy for V^ 1Hg and 
Vj^ lHi , even after taking the different assumptions noted above into account. Our efforts to 
reconcile the differences have been unsuccessful. 

Gunaydin et al. have performed first-principles calculations for neutral Vaih 3 and Vn & h 
in a 192- atom NaAlH4 supercell. js] Their calculated formation energies are in reasonably 
good agreement with our present values, for comparable choices of chemical potentials. They 
obtained free-energy barriers for diffusion using an umbrella-sampling technique, resulting 
in values of 0.12 eV for Vaih 3 and 0.26 eV for Vn & h- We feel that these computed free-energy 
differences are not representative of actual diffusion activation energies, which should only 
include enthalpy differences and do not depend on entropy terms in the free energy. This 
explains the difference between the values reported by Gunaydin et al. and our values listed 
in Table [H 

We also note that Gunaydin et al. 9| seem to have scaled their calculated formation en- 
ergies and migration barriers by a factor 2/3, presumably to convert from units of "kJ/mol" 
to "kJ/mol H 2 " because in Eq. (3) of their paper 3/2 H 2 is being produced. Such scaling is 
unwarranted and unjustified, in our opinion. The formation energies and migration barriers 
enter into activation energies, which do not scale with the size of the system. Put differ- 
ently: producing a larger amount of H2 requires a given process to happen a larger number 
of times, but it does not increase the activation energy of that process. In the absence of this 
scaling factor, i.e., looking at their directly calculated numbers, the conclusions of Gunaydin 
et al. are not supported by their calculations. 

B. Migration barriers 

As reported in Sec IHIt migration energies were explicitly calculated in the cases of V# , 
Vft, H~, VJ^ a , Na^~, V^J H4 , and (H 2 )j. Explicit calculations for migration paths proved 
very cumbersome in other cases, specifically H+, V^an, Vaih 3) Vav an( ^ ^at The rea son 
for the complications is that these defects can be thought of as complexes between more 
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elementary native point defects, as discussed in Sec. IIHI In this view, however, one can also 
estimate a lower bound for the migration energies of these complexes by taking the higher 
of the migration barriers of the constituent defects. For instance, can be considered 
as a complex of and (H 2 )j, suggesting that its migration barrier will be at least 0.26 
eV, the value for V^. This estimate should be a lower limit for the migration barrier of 
the complex, since it assumes that during the migration process the complex remains fully 
bound. If during migration (partial) breaking of the complex occurs, then an energy cost 
needs to be paid that would increase the value of the barrier. Values for other defects 
that can be regarded as complexes were obtained in similar fashion: the aluminum-related 
vacancies can be regarded as combinations of V^J H , (H 2 )j, and H^". Vn & h is regarded as a 
combination of K7 and Kt. 



C. Reaction mechanisms and activation energies 



A range of apparent activation energies for desorption in transition-metal-doped NaAlfL 
has been reported, with values as low as 0.8 eV. This value is ~0.4 eV lower 

than in undoped material. 3fl The desorption process has been suggested to arise from the 

nn 

diffusion of defects. |7|, |8J 

In order for hydrogen to desorb from this material, it is also necessary for aluminum 
and/or sodium to diffuse through the solid. The reaction described by Eq. ([I]) can be in- 
terpreted as the material NaAlH4=(NaH)(Alfi3) decomposing into Na3AlH6=(NaH)3(Alfi3) 
and AlH3=(NaH)o(AlH3), with the latter divided into two separate phases - Al and H2. 
This leads us to be interested primarily in the diffusion of species corresponding to fluxes of 
either (NaH) or (A1H 3 ) through the material. 

In Sec. IHI Al we explained that for purposes of presenting our formation energy results, 
we chose the chemical potentials of Al, Na, and H by assuming equilibrium with NaAlH 4 , 
NasAlHg, and Al. Based on the calculated enthalpies of formation, and keeping in mind 
that the chemical potentials are referenced to the standard state of the elements, this leads 
to values /iNa=— 0.38 eV, /iai=0, and /ih=— 0.12 eV. The resulting formation energies were 
presented in Figs. [2j [H and El and in column (1) of Table [B Of course, other scenarios are 
possible. For instance, we can assume equilibrium with NaAlH 4 , Na3AlH 6 , and H 2 , and find 
/iNa= — 0.48 eV, /iAi=— 0.34 eV, and /ih=0. The resulting formation energies are listed in 
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column (2) of Table [U with the additional assumption that again the Fermi level is fixed 
at the value where the formation energies of V^J H and are equal. As a third scenario, 
we can assume equilibrium with NaAlH 4 , Al, and H 2 , finding /i Na =— 0.82 eV, /iAi=0, and 
/iH = 0. The resulting values are listed in column (3) of Table [B 

We are now in a position to examine various possible mechanisms for defect-assisted 
diffusion and decomposition. The reaction described by Eq. (JT|) implies formation of Na3 A1H\6 
and Al during the decomposition of NaAlH 4 . Since the formation energies depend on atomic 
and electronic chemical potentials, specific choices must be made that approximate the 
experimental conditions as closely as possible. Our assumption is that these conditions do 
not differ significantly from equilibrium. In the case of charged defects, local and global 
charge neutrality needs to be maintained. 

1. Neutral defects 

One possible mechanism involves diffusion of neutral defects such as Vaih 3 or VneH- The 
formation energy of these defects is independent of Fermi level. The total activation energy 
would be the sum of formation and migration energies. Vaih 3 and VNa cannot form in 
the interior of the material, since this would require simultaneous formation of Alj or Na^, 
and the total process would be much too costly. Al- and Na-related defects will therefore 
necessarily be formed at an interface, where they will be directly involved in the formation of 
NaaAlHg and/or Al. Vaih 3 most readily forms at an interface where NaAlH 4 , Al and H 2 are 
in equilibrium (as also discussed in Ref. 0). Table [J shows that under those conditions the 
activation energy for Vaih 3 diffusion is given by 1.13+0.46=1.59 eV. Creation of VNaH would 
happen at an interface between NaAlH 4 and NasAlHg, resulting in an activation energy of 
1.25+0.41 = 1.66 eV. These values are actually lower bounds, since as explained in Sec. IIVBI 
the migration barriers used here are only estimates. However, even these lower bounds are 
already higher than the observed activation energies for undoped alanate. We therefore 
conclude it is unlikely that neutral Vaih 3 or Vn & h entities would play a dominant role. 

As pointed out above, these defects are actually not elementary defects in their own 
right. Vaih 3 should be regarded as a complex consisting of V^ H4 and H,^, while V^aH is 
a combination of and V^. Indeed, based on the structure of NaAlH 4 it is clear that 
V^ih and V^~ a should be regarded as the elementary Al- and Na-transporting entities in this 
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material; note, e.g., the persistence of Al-H vibrational modes in NaAlH4 up to the melting 
temperature. [23] 



2. Other Al- and Na-related defects 

Table [I] shows that it is too costly for defects containing only Al to diffuse. The estimated 
activation energies, both for Vai and for Alj are simply too high. For Na-related defects, V^ a 
is clearly important and will be discussed in more detail below. Na^~, on the other hand, is 
too costly to form. 

3. Role of hydrogen- related defects 

In Sec. IIV C II we noted that V^ lHi and VJ^~ a should be regarded as the elementary de- 
fects responsible for Al and Na transport. Note that we do not necessarily need both to 
explain the decomposition of NaAlH^; either defect by itself would suffice to explain forma- 
tion of NasAIHe and Al. The accompanying defects (H~ and V^) discussed in Sec. IIV C II 
are hydrogen-related defects that provide local charge neutrality and additional hydrogen 
transport; however, there is no pressing reason why, out of possible hydrogen-related defects, 
Hr should be the only one associated with V^ m4 , or the only one associated with V^~ a . 
The formation energy of the hydrogen-related defects is sensitive to the choice of chemical 
potentials and the position of the Fermi level, and different conditions may favor differ- 
ent defects. Table U shows that among the positively charged defects, is systematically 
lower in energy than H+, and among the negatively charged defects, H~ is lower in energy 
than Vjf. However, different Fermi- level conditions will favor either positively or negatively 
charged defects. 

It is then important to realize that, from the point of view of transporting hydrogen, 
is equivalent to V^. Indeed, moving a hydrogen interstitial from left to right is equivalent to 
moving a vacancy from right to left. And given the opposite charge on these defects, local 
electric fields would indeed tend to push these defects in opposite directions. 

We look upon H 4 ~ and as charge- carrying defects that will provide the necessary local 
charge neutrality. In the absence of electronic charge carriers (electrons or holes), it is these 
highly mobile hydrogen-related defects that will provide local charge balance in an insulator 
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such as NaAlH4. 

Unlike the Na- or Al-related defects, the hydrogen-related defects can readily form in the 
interior of the material. H~ and V^" constitute a Frenkel pair that can be formed simply 
by moving a substitutional H atom to an interstitial position. This contrasts sharply with 
the Al- or Na-related defects, which cannot form in the bulk but only at an interface, as 
discussed in Sec. II VC II The hydrogen-related defects are thus unique in their ability to form 
within the bulk, and given their modest formation energies, we expect a finite concentration 
of such pairs to always be present. 

We suggest here that it is the formation of these hydrogen defect Frenkel pairs that may 
be the rate limiting step for the desorption and decomposition process. Our calculated 
formation energy for the Hj-Vfj pair is 1.10 eV, corresponding to the configuration where 
the pair creates an AIH4-H-AIH3-H-AIH3 (or AI3H12) complex. Note that this is lower than 
the sum of the formation energies of H~ and listed in Table (TJ the reason being that the 
proximity of the defects leads to a finite binding energy, which we calculate to be equal to 
0.46 eV. This calculated formation energy would yield an estimate for the activation energy 
of formation of this pair of 1.10+0.16=1.26 eV, where we have added the migration barrier 
of to the formation energy as an estimate for the barrier that needs to be overcome. 
This activation energy is consistent with the experimentally observed desorption activation 
energies in undoped NaAlH4. 

This suggested mechanism also immediately provides an explanation for the observed 
effect of Ti (or other transition metals), following the arguments in Ref. 2). Ti acts as an 
electrically active dopant that shifts the Fermi level. Shifts of the Fermi level will always 
lead to a lowering of the formation energy of the hydrogen-related defects. Q For purposes 
of illustration, consider an upward shift in the Fermi level, such as would be introduced by 
Ti. 2! Figure [2] shows that if the magnitude of this shift exceeds a few 0.1 eV, it becomes 
energetically more favorable for the vacancy to be present in the negative charge state rather 
than the positive charge state. The H~-V^~ Frenkel pair will then be lower in energy than 
the H^-Vg" pair. To illustrate the effect of a shift in Fermi level, consider the example 



given in Ref. 
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I an upward shift of 0.44 eV in the Fermi level leads to a decrease in the 
formation energy of by 0.44 eV and an increase in its concentration [Eq. (j3j)] by 6 orders 
of magnitude at 100°C We suggest that this energy gain, in the presence of an electrically 
active impurity that shifts the Fermi level, explains the lowering in activation energy for 
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decomposition and desorption. Note that the components of the -V K pair repel each 
other, which eliminates any further energy cost in separating them into isolated defects. 

V. CONCLUSIONS 

We have studied native defects in NaAlH4 based on first-principles density functional 
calculations. Despite the wide variety of possible native defects in this material, simplifying 
principles can be applied to qualitatively understand the relevant defects. Many of the 
defects can be described as complexes containing the following elementary defects in the 
system: V^ H4 , V^ a , V^, H~, and (H 2 )j. This does not eliminate the need for performing 
explicit calculations, since the interactions between constituent defects significantly lower 
the energy of the complex, e.g., through relaxation of local strains. We suggest that this view 
of native defects will also be useful when studying other alkali or alkaline earth aluminum 
hydrides. 

Our results should be of use in interpreting experimental information on desorption from 
NaAlH 4 and the effects of dopants such as transition metal impurities on the kinetics. 
Charged (rather than neutral) defects play the dominant role in diffusion and decomposi- 
tion. We have proposed specific mechanisms for this process; in particular, we suggest that 
hydrogen-related Frenkel pairs provide the charge carriers that aid in diffusion of Al- or Na- 
related defects. The dependence of their formation energy and, therefore, their concentration 
on Fermi level may explain the observed effects of transition-metal dopants. Moreover, one 
could engineer the hydrogen desorption kinetics by adding electrically active impurities to 
shift the Fermi level for tuning the concentrations of the defects relevant for hydrogen des- 
orption. The impurities suitable for that purpose do not necessarily need to be transition 
metals. 
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TABLE I: Calculated formation energies E f and migration energies E m for selected defects in 
NaAlH4. For charged defects, the formation energies are taken at the Fermi- level position where the 
formation energies of V^ m ^ and are equal. (1), (2), and (3) refer to different choices of chemical 
potentials. (1) corresponds to equilibrium with NaAltLi, NasAlHg, and Al (/iN a =— 0.38 eV, /iAi=0, 
^=-0.12 eV; e F =3.26 eV); (2) to equilibrium with NaAlH 4 , Na 3 AlH 6 , and H 2 O Na =-0.48 eV, 
^Ai =— 0.34 eV, /Uh=0; e^=3.15 eV); and (3) to equilibrium with NaAlELi, Al, and H2 (//Na=— 0.82 
eV, ^ai=0, and nu=0; ep=2.81 eV). Conditions (1) correspond to the case depicted in Figs. El [U 
and Migration energies denoted by an asterisk (*) are lower bounds estimated by considering 
the defect as a complex and taking the higher of the migration energies of the constituents; see 
text. 



Defect 




Ef (eV) 




E m (eV) 




(1) 


(2) 


(3) 




1 11 


0.90 


0.91 


0.57 


0.26 


Vh 


1.00 


1.23 


1.57 


0.92 


H+ 


1.27 


1.04 


0.70 


0.26* 


HT 


0.66 


0.65 


0.99 


0.16 


(H 2 )i 


1.07 


0.84 


0.84 


0.25 


v& 


1.12 


0.67 


0.67 


0.46* 




1.13 


0.90 


1.58 


0.46* 


V A\ 


1.67 


1.66 


3.02 


0.46* 


V+ 
t/ AlH 4 


0.49 


0.50 


0.50 


0.46 


V° 
AIH3 


0.78 


0.79 


1.13 


0.46* 


^A1H 2 


1.07 


0.96 


1.30 


0.46* 


Alf 


4.18 


4.19 


2.83 






0.49 


0.50 


0.50 


0.41 


^NaH 


1.25 


1.27 


0.93 


0.41* 


Na+ 


1.44 


1.43 


1.43 


0.48 



18 




FIG. 1: (Color online) The body-centered tetragonal structure of NaAlELi. Sodium is represented 
by large grey spheres, aluminum (at center of tetrahedra) is small dark blue, and hydrogen is small 
red. A [001] plane is drawn in grey in order to help illustrate depth. 
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Fermi level (eV) 

FIG. 2: Formation energy of hydrogen-related native defects in NaAlH4, plotted as a function of 
Fermi level with respect to the valence-band maximum. Atomic chemical potentials were chosen 
to reflect equilibrium with NasAlHg, NaAlH4, and Al. 
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FIG. 3: (Color online) The structures of (a) Hr, (b) V+ , (c) V R , (d) (H 2 );, and (e) H+. In the 
case of H~ and an isosurface of the highest occupied orbital is shown, in order to illustrate the 
defect-character of this orbital. The electron density at the isosurfaces is 0.06 electrons/A 3 . All 
of the structural diagrams are of sodium and aluminum atoms in a [010] plane (and coordinated 
hydrogens), except for the case of (H^j in panel (e) where a [001] plane is shown. 
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2 2,5 3 3.5 4 



Fermi level (eV) 

FIG. 4: Formation energy of aluminum-related native defects in NaAlH4, plotted as a function of 
Fermi level with respect to the valence-band maximum. Atomic chemical potentials were chosen 
to reflect equilibrium with NasAlHg, NaAlELi, and Al. 
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FIG. 5: (Color online) The structures of (a) V^J H4 , (b) V^, (c) V^ - , (d) V^j, and (e) Vaih 3 - The 
structural diagrams in panels (a), (b), and (e) are of sodium and aluminum atoms in a [010] plane 
(and coordinated hydrogens). For the cases of and (c) and V^j (d) , we show the atoms from 
two [010] planes in order to illustrate the differences upon reduction of the aluminum vacancy. 
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FIG. 6: Energy versus reaction coordinate along the migration path for V^J H4 (equivalently, A1HJ) 
as modeled with the NEB method. The local lattice structures of the initial, saddle-point, and 
final configurations are shown as (I), (II) and (III), respectively. 
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Fermi level (eV) 

FIG. 7: Formation energy of sodium-related native defects in NaAlH4, plotted as a function of 
Fermi level with respect to the valence-band maximum. Atomic chemical potentials were chosen 
to reflect equilibrium with NasAlHg, NaAlH4, and Al. 
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FIG. 8: (Color online) The structures of (a) and (b) V^aH- The structural diagrams are of 
sodium and aluminum atoms in several [010] planes (and coordinated hydrogens). 
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